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We present a construction for improving numerical cubature formulas with equal weights and a convolution 
structure, in particular equal-weight product formulas, using linear error-correcting codes. The construction 
is most effective in low degree with extended BCH codes. Using it, we obtain several sequences of explicit, 
positive, interior cubature formulas with good asymptotics for each fixed degree t as the dimension n — > °°, 
Using a special quadrature formula for the interval 1 13], we obtain an equal-weight f-cubature formula on the 
n-cube with 0(«L'/ 2 J ) points, which is within a constant of the Stroud lower bound. We also obtain f-cubature 
formulas on the ^-sphere, n-ball, and Gaussian R" with 0(n'~ 2 ) points when t is odd. When /i is spherically 
symmetric and t = 5, we obtain 0(n 2 ) points. For each t > 4, we also obtain explicit, positive, interior formulas 
for the n-simplex with 0(n'~ l ) points; for t = 3, we obtain 0(n) points. These constructions asymptotically 
improve the non-constructive Tchakaloff bound. 

Some related results were recently found independently by Victoir l2lll . who also noted that the basic con- 
struction more directly uses orthogonal arrays. 



1. GENERAL RESULTS 

Let ji be a normalized measure on M." with finite moments. 
A cubature formula of degree t, or t -cubature formula, for jx 
is a set of points F = {p a } CU° and a weight function p a i— > 
w fl £ K such that 

r N 

P(x)dii=P(F)=Y< w « P (P«) 

J a=\ 

for polynomials P of degree at most t . (If n = 1, then F is also 
called a quadrature formula.) The formula F is equal-weight 
if the w a are all equal; positive if w u > for all a; and other- 
wise it is negative. Let X be the support of fi. The formula F 
is interior if every point p a is in the interior of X; it is bound- 
ary if every p a is in X and some p u 6 dX; and otherwise it is 
exterior. These properties of cubature formulas are often ab- 
breviated. E.g., PI means positive and interior and EB means 
equal-weight and boundary. (Exterior formulas are denoted 
"O," for outside.) An equal-weight formula is abbreviated "E" 
and is also called a (geometric) t-design or a Chebyshev-type 
formula. 

The main use of a cubature formula is to numerically inte- 
grate a function / which is approximately a polynomial. In 
this application, formulas with many points or non-explicit 
points are impractical, exterior formulas are ill-founded if / is 
only defined on X, and formulas with large negative weights 
are ill-conditioned on the class of continuous functions 12(1 
Ch. 1]. Thus PI formulas with few points are the best kind. 

By Tchakaloff's theorem [2Cj p. 61], every measure fi on 
R" has a PI f-cubature formula with at most \V) points, the 
same as the dimension of the vector space of relevant polyno- 
mials, M[3c]< f . (If dX has non-zero measure, it may only be 
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a PB formula.) Tchakaloff's theorem has a short proof, but 
it is computationally non-constructive. Many known formu- 
las with n small, or with n large and t < 2, are better than 
the Tchakaloff bound l5l l20ll . But if n is large, t > 3, and jj, is 
reasonably natural, most explicit formulas in the existing liter- 
ature are either negative, exterior, or have exponentially many 
points. 

In this article we present a new method to thin equal-weight 
cubature formulas with a convolution structure, in particular 
product formulas for product measures. (By thinning a for- 
mula, we mean removing some of its points without reducing 
its cubature degree.) The thinned formulas are efficient in high 
dimensions and low degree. The method also applies to some 
non-product measures that are related to product measures, in 
particular spheres and simplices with uniform measure. Vic- 
toir 12 ill independently obtained the basic construction when 
q = 2, together with some other generalizations not consid- 
ered by this author. However, many of our asymptotic bounds 
and derived constructions are new. 

If F and G are two cubature formulas, we define their con- 
volution F * G to be their sum as sets, F + G. The weight w a 
of p u in F * G is given by a product rule: 

Pa=Pb+Pc 

Pi,eF,p c beG 

Convolution of cubature formulas is related to convolution of 
measures in two ways: First, it is convolution of measures if 
cubature formulas are interpreted as atomic measures. Sec- 
ond, if F is a f-cubature formula for fL and G is a f-cubature 
formula for v, then F * G is a f-cubature formula for fl * V. In 
particular, product formulas and product measures are convo- 
lutions in independent directions. 

We also recall some basic facts from coding theory. For 
each prime power q, there is a unique finite field ¥ q with q el- 
ements. A linear error-correcting code of length £, dimension 
k, and distance t over F 9 is a fe-dimensional vector subspace 
of W„ such that each non-zero vector has at least t non-zero 
coordinates. It is also called an [l,k,t] a code. A code C is a 
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zero-sum code if the coordinates of every a G C sum to 0. 

Theorem 1.1. Let t, n, and £ be positive integers, let q be a 
prime power, and let jl be a measure on W. For each 1 < i < 
£, let Fj be an equal-weight formula with q elements such that 
the convolution 

F = F\ * F% * . . . * Fi 

is a t-cubature formula for /X. Then an [£,k,t + \] q code C 
yields a thinning G C F with q f ~ k points. In addition, if each 
Fi is centrally symmetric, t is odd, and either q is odd or C is 
a zero-sum code, then C need only be an [£,k,t]g code. 

Theorem 1 1.1 1 can be strengthened further using the notion 
of an orthogonal array . Linear error-correcting codes are 
dual to linear orthogonal arrays, and the proof actually uses 
orthogonal arrays rather than codes. In some cases non-linear 
orthogonal arrays are slightly better than linear ones. See Sec- 
tions|§]and|4] 

The most effective case of Theorem ll.ll is in the asymptotic 
limit n — > °° with t and q fixed. Recall that a function /(«) is 
quasilinear if f(n) = O((logn) a n) for some a. Quasilinearity 
is also written /(«) = 0(n). Say that a family {F} of cubature 
formulas is quasilinear (abbreviated "QL") if the points and 
weights of each F can be generated in quasilinear time in the 
length of the output. 

Theorem 1.2. Assume all variables as in Theorem U.h Then 
G can have 0(£ a ) points (with the constant depending only 
on q), where 

, r ^ — 1 i 
a = r- 1 - . 

q 

If each G is centrally symmetric and t is odd, then 

„ i f - 2 i 

a = t-2- . 

q 

Moreover, G is quasilinear as I — > °° assuming precomputa- 
tion of each Fj. 

If jj. is an ra-fold product with m n, then we can take £ °< n 
in Theorem 1 1.21 so that 0(£ a ) = 0(n a ). In comparison, the 
Tchakaloff upper bound is 0(n') points, or 0(n'~ l ) when t 
is odd and fj, is centrally symmetric (Section^}. Thus Theo- 
rem ll.2l is asymptotically better than Tchakaloff 's theorem for 
all such product measures. Tchakaloff 's theorem also does not 
guarantee equal weights. 

Another comparison is with the cardinality of exact de- 
termination. A f-cubature formula F is overdetermined, un- 
derdetermined, or exactly determined if the parameters of its 
points provide fewer, more, or the same number of degrees 
of freedom, respectively, as the constraints imposed by inte- 
grating all polynomials of degree t. The cardinality of exact 
determination is ©(«' _1 ) for general jl and 0(n' _1 ) when t is 
odd and jl is centrally symmetric. Thus for product measures, 
the formulas in Theorem 1 1.21 are asymptotically exactly de- 
termined (up to a constant factor that depends on t ) when q is 



large. But when q<t— 1, or q<t — 2 in the odd and centrally 
symmetric case, they are asymptotically overdetermined. 

A third comparison is with the Stroud lower bound: Any f- 
cubature formula in n dimensions, not necessarily interior or 
positive, requires n(nL'/ 2 J) points. Theorem l 1 .21 achieves the 
Stroud bound (up to a constant factor) when q = 2. 

A final comparison is with an interesting thinning construc- 
tion of Novak and Ritter for products of quadrature formulas 
fl^l . (It is similar to an earlier construction due to Grundmann 
and Moller for the n-simplex |7].) They produce f-cubature 
formulas with <9(«L'/ 2 J ) points, which is within a constant fac- 
tor of the Stroud bound and better than Theorem 11.21 when 
q > 2. Crucially, their formulas are not positive, although 
they can be made interior. They also require that the factors 
of jl be 1 -dimensional. The Novak-Ritter construction does 
generalize to convolutions, as long as each factor formula has 
collinear points. 

Theorem 1 1.21 can be used to construct interesting cubature 
formulas for several infinite sequences of regions and mea- 
sures considered by Stroud 1 20, Ch 7,8]: 

Theorem 1.3. For any t: 

1. The n-cube C n with uniform measure has a QLEI t- 

cubature formula with 0(«L'/ 2 J ) points. 

2. The cubical shell C n — rC„ has a QLEI t-cubature formula 

with 0(nL'/2J+ 1 ) points. 

For any odd t > 3: 

1. R" with Gaussian weight function has a QLEI t-cubature 

formula with 0(n'~ 2 ) points. 

2. Any spherically symmetric measure on W has a QLPI t- 

cubature formula with 0(n'~ 2 ) points. This includes the 
n-ball B n , the spherical shell B„ — rB„, and the (n— 1)- 
sphere S"~ 1 with uniform measure; and R" with radial 
exponential weight function exp(— | \x\ I2). 

For any t > 2: 

1. The n-simplex A„ has a QLPI t-cubature formula with 

(9(n f_1 ) points. 

2. The n-cross-polytope C* with uniform measure has a QLPI 

t-cubature formula with <9(« L 3r / 2 J — 1 ) points. 

All cases of Theorem l 1 . 3l other than the cross polytope C* 
improve the Tchakaloff bound. On the other hand, the con- 
struction for the cube C n matches the Stroud bound up to a 
constant factor. We admit that this f-dependent factor is very 
generous when t is large: For each t = 2s + 1, it approaches 
2 • s s ■ s\ as n — > °° in the favorable case n — 2 m . By contrast, 
the Novak-Ritter formulas use only 2 s more points than the 
Stroud bound as n — > °°. 

Theorem 11.31 partially solves a problem of Stroud J20L p. 
18]: Are there PI 5-cubature formulas for C„, B„, or A„ 
with 0(n 2 ) or 0(n 3 ) points? Theorem ll.3l provide s QLPI 5- 
cubature formulas with 0(n 2 ) points for C„, <9(n 3 ) points for 
B n , and 0(n 4 ) points for A„. In Section|3j we will establish 
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a special QLPI 5-cubature formula for B„ with 0(n 2 ) points 
and QLPB and QLPI 3-cubature formulas for A„ with 0(n) 
points. Thus the only remaining case of Stroud's question is 
the ^-simplex in degree 4 or 5. 

Remark. The formula in Theorem ll.3l for S"~ 1 is technically a 
QLPB formula if we take the definition of boundary in general 
topology. However, we take boundary in the sense of geomet- 
ric topology, so that Theorem ll.3l is correct as stated. 
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2. PROOFS 

Proof of Theorem ]! .1\ First, identify an affmely independent 
set of q points in W~ l with the finite field ¥ q . For each 1 < 
i < A choose a linear map %i : W 1 - 
and define K : R^ -1 ^' — > R" to be their direct sum: 



that sends ¥ q to Ft, 



71 = 7Ti © 7T2 ( 



) 7Z{. 



Because F\ *Fz*- ■ ■ *F[ is a t -design for the measure fi on R", 
the identity 

J q P^n 

holds for any polynomial P of degree at most t on R". Now 
suppose that we thin the set F = Tr(F^) to a set G — 7c(A) 
for some set A C Wi. Since % is linear, if we want G to be a 
f -cubature formula for \i as F is, it suffices that 



q pe¥' a 11 peA 



(1) 



for any polynomial P on MS q l ^ of degree at most t. If P is 
a monomial, then as a function on ¥ q it depends on at most 
t coordinates. Conversely, any function on F^ it depends on 
at most t coordinates is realized by a polynomial of degree at 
most t. It follows that equation Q is equivalent to the sta- 
tistical property that the projection of A onto any t of the / 
coordinates of ¥„ is constant-to-1. Such a set A is called an 
orthogonal array of strength t. 

If C is an [£,k,t + l] q code, then the dual space C* (in the 
sense of linear algebra over ¥ q ) is a linear orthogonal array of 
strength t. Since C has dimension k, C* has dimension £ — k 
and therefore has gr* points. Thus we can let G = 7t(C*). 

The refinement when t is odd and each Fi is centrally sym- 
metric is as follows. If q is odd, we replace W~ l by R^ -1 )/ 2 , 
and we position ¥ q as a centrally symmetric set that does not 
lie in a hyperplane. (In other words, the points of ¥ q are the 



vertices of an affmely regular cross polytope, plus the origin.) 
We further demand that negation in ¥ q coincides with nega- 
tion in R^ -1 )/ 2 . Then any centrally symmetric subset A C F^ 

is centrally symmetric in MS q ~ 1 ^' 2 . In this case both sides of 
([0 vanish when P is an odd polynomial. Thus A need only 
be an orthogonal array of strength t — 1. In particular, this is 
so if A = C* , because C* is a vector space over ¥ q and vector 
spaces are centrally symmetric sets. 

Finally if t is odd, q is even, and C is a zero-sum code, then 
C* contains the vector (1, 1, . . . , 1) and is therefore invariant 
under addition by this vector. In this case we replace W 1 ^ 1 in 
the general construction by R 9 / 2 and we realize ¥ q as a cen- 
trally symmetric set (the vertices of a regular cross polytope). 
We further demand that adding 1 in ¥ q coincides with nega- 
tion in R ? / 2 . Then once again C* is centrally symmetric and 
need only be an orthogonal array of strength t — 1 . □ 



The following lemma establishes Theorem 1 1.21 as a corol- 
lary of Theorem ll.il 

Lemma 2.1. Let q be a prime power, let m,t £ Z>o, and let 

o=*-i- l— j. 
i 

Then there is a [q m ,k,u] q zero-sum code C with 
u>t + l k>q m - ma - 1 . 

The code in Lemma [2.11 is called an (extended, narrow- 
sense) BCH code H HUE Q- We will use the duals of 
BCH codes to thin cubature formulas. As it happens, the dual 
of a BCH code of this type is another BCH of the same type. 

Proof. It is easier to define the dual code C* and show that it 
is an orthogonal array. Since it is a linear space, it suffices to 
show that every coordinate projection Tij : C* — > F^ with |/| < t 
is onto. There is an important F^-linear function 



Tr 9 : ¥qm 



called the trace. (It is analogous to the taking the real part of 
a complex number.) First, we interpret F c ; as the space of 
all functions from ¥ q m to F^. We define C* as the set of all 
functions 



f(x)=TT q (P(x)), 



where P is a polynomial of degree at most t — 1 . If / C ¥ q m 
and |/| < t , the polynomial P can achieve any desired values 
on / by Lagrange interpolation. Thus the distance of C is at 
least t + 1 . 

The space of polynomials of degree t — 1 on ¥ q m has F^™- 
dimension t, and therefore F^-dimension mt. But taking the 
trace reduces the dimension in two ways. To give an explicit 
example, suppose that q — 2, t = 3, and m is arbitrary. Then 
C* is the set of all 

f(x) =Tv 2 (ax 2 + bx + c). 
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The apparent dimension of C* is 3m. But / only depends on 
the trace of c, so c contributes 1 rather than m to the dimension 
of C* . Moreover, Tr2(fct) = Tx2{b 2 x 2 ), so the linear term can 
be removed from /, with the conclusion that 

dimC* < m + 1 . 

In general, the constant term of P contributes 1 to the dimen- 
sion and the other t — 1 terms contribute m each, except that 
I —J terms are superfluous by the Frobenius automorphism 
x i— > x q . Thus 

dimC < ma + l, 

as desired. 

Since constants are polynomials of degree 0, C* contains 
constant vectors. Therefore C is a zero-sum code. □ 

On the face of it, Lemma l2~Tl onIv establishes Theorem ll.2l 
when £ = q m . If q"'~ l < £ < q' n , we can project a BCH code 
from ¥ H q to F^. This preserves the (9(f ) bound at the ex- 
pense of worsening the constant factor. If £ is not much more 
than q" 1 ^ 1 , we can slightly improve the projected code with a 
projection that annihilates up to a — 1 independent vectors in 
C*. (See Theorem O.ll for an example.) 

Remark. The inequalities for u and k in Lemma l2~Tl become 
sharp as m — > °°. 

Proof of Theorem ]! .31 The simplest case to consider is with 
uniform measure and W with Gaussian weight function. This 
fits Theorem ll .21 with £ = n, provided that for each f, we find 
an EI f -quadrature formula with Gaussian weight and with q 
points for some prime power q. Since there is no bound on q, 
the Seymour-Zaslavsky theorem fl^l establishes that such for- 
mulas exist. One explicit method begins with the PI Gaussian 
(2f + l)-quadrature formula with t + 1 points. Viewed as a t- 
quadrature formula, the t + 1 points can be freely perturbed. 
In particular they can be perturbed so that the weights become 
multiples of 1/q for some large prime power q. The pertur- 
bation can be chosen to retain central symmetry. On the other 
hand, since q is large, [(t — 2)/q\ —0. Thus Theorem 1 1.21 
produces formulas with 0(n l ~ 2 ) points. 

We will need the same construction for the orthant R" with 
exponential weight function exp(— This measure does 

not have central symmetry, and the end result is formulas with 
0(n ) points, again with £ — n. 

The next simplest case is the n-simplex A„. Recall that A„ 
has barycentric coordinates 

Xq+X\ + ... +x„ = 1 

which realize it as a subset of the orthant R" q 1 . If P(x) is a 
polynomial of degree t on A„, then it can be homogenized: it 
can be expressed as a homogeneous polynomial of degree t 
by attaching a factor of (L,x,) f S to each term of degree s. In 
this case 

/ P{x)dx = 1 j P{x) exp(- 11*11! )dx. 
Ja„ {n+ty.Jw!* 1 



Therefore we can project any non-exterior cubature formula 
for R"^ 1 radially onto A„ without loss of degree, although 
the weights change. (If the origin happens to be a cubature 
point, discard it.) In particular, we can project the cubature 
formulas provided by Theorem 1 1.21 as explained previously. 
The formulas still have 0{n'~ l ) points, although the weights 
are no longer equal. 

The same argument works for the sphere S"~ 1 C R" for cen- 
trally symmetric formulas. Every polynomial P on S"~ l can 
be expressed as P$ + Pa, where P$ is centrally symmetric and 
Pa is centrally antisymmetric. The integral of Pa vanishes, 
as does its sum with respect to any centrally symmetric for- 
mula. Meanwhile every term of P$ has even degree, so it can 
be expressed as a homogeneous polynomial on R" using the 
equation 

x} + ...+x 2 „ = 1 
for the unit sphere. Then 

/ P{x)d& - 2 / P(x) exp(- 1 \x\ \f)dx, 

where Q. is usual surface volume on S"~ 1 . Again, any centrally 
symmetric cubature formula can be radially projected and the 
weights adjusted. 

Formulas for the ball B n and the spherical shell B„ — rB„ 
can be derived from formulas for the sphere S"~ l using radial 
separation of variables I20l Th 2.8]. The result is a product 
formula where the radial factor can be Gaussian quadrature. 
The number of points in this factor does not increase with 
dimension. 

The cross -poly tope C* is the union of 2" simplices. Thus 
we can obtain formulas for C* by repeating formulas for A„. 
In degree t , we do not need all 2" copies; instead we can repeat 
it in the pattern of the BCH code over F2 defined by polyno- 
mials of degree t — 1 over V%m. Such a code has 0(nL ; / 2 J) 
vectors and the formula for A„ has 0(n'~ 1 ) points, so the total 
is 0(n L3r /2J-i) points ^ 

The n-cube C„ — [— 1 , 1]" is in some ways the most interest- 
ing case. Like the Gaussian case, it is a straight application of 
Theorem ll. 21 using an equal-weight quadrature formula. But 
in this case we will carefully choose the quadrature formula 
on [—1,1] to itself be a convolution of s = [t/2\ formulas 
with two points. For example, the Chebyshev 5 -quadrature 
formula has points at 

I5+V5, I5-V5 

This is evidently a convolution, as is any centrally symmetric, 
equal-weight formula with 4 points. Elsewhere I13il we show 
that the 2 s points 

±Z1± Z2 ±...±Z. 

form a Chebyshev-type (2s + 1) -quadrature formula for 
[—1,1] with constant weight if and only if the z,'s are the roots 
of the polynomial 

v x s - 1 X s - 2 (-1V 
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We also show that all roots of Q are real and that the resulting 
quadrature formula is interior. The n-fold product power of 
this formula is thus a convolution of sn pairs of points, so we 
can apply Theorem l 1 .21 with i = sn and q = 2. 

Finally the 0(n^'^) formula for the «-cube C„ yields a 
0(«L r / 2 J+i ) formula for the cube surface dC„ just by repeating 
the formula for C„_i on each facet of C„. Then radial sepa- 
ration of variables produces a product formula for the cubical 
shell C„ — rC n which also has (9(«l-'/ 2 J+ 1 ) points. □ 



Kerdock OA(2 4m ,2 2m ,2,5), and then projected down to n di- 
mensions. 

If 2 2m <n< 2 2m + 2"\ then the [2 2m+l ,4m + 3,6*] 2 BCH 
code can be reduced by half by carefully choosing the projec- 
tion. The code has a vector of weight 2 2m — 2™, so when n is 
only slightly larger than 2 2m , we can choose a projection that 
annihilates this vector. 

In each of the three cases, the result is a formula with 2 k + 
2n points. □ 



3. SPECIAL CONSTRUCTIONS AND EXAMPLES 

In this section we will consider some examples and spe- 
cial constructions with concern for constant factors. For this 
purpose, we spell out more precisely the notion of an orthog- 
onal array. Let A be a finite set. If a subset X C A" has the 
property that its projection X — > A ! is a constant-to-1 map for 
every |/| < t, then T is an orthogonal array of strength t , or an 
OA(|r|,n,|A|,f) 191 If A = F 9 and X = C* is the dual of an 
[n,k,t]q code, thenX is an OA(q"~ k ,n,q,t — 1). We will also 
say that X is an [n, n — k, t*] q to refer to its linear structure and 
indicate its dual distance. 

If \S\ — q is a prime power and t is fixed, then BCH codes 
are the best presently known F^-linear orthogonal arrays in the 
limit n — > °°. But a few non-linear arrays are slightly better. 

A Hadamard matrix of order n is an n x n matrix with en- 
tries ±1 and with orthogonal rows (and therefore orthogonal 
columns as well). It is easy to show that a Hadamard matrix is 
equivalent to an OA (2n , n , 2 , 3 ) . A [2 m , m + 1 , 4*] 2 BCH code, 
which is also called a first-order Reed-Muller code, yields a 
Hadamard matrix of order 2 m . But there are also Hadamard 
matrices for other values of n, for example when 4\n and n — 1 
is prime. The Hadamard conjecture asserts that there is a 
Hadamard matrix of every order n divisible by 4. 

For any even m > 4, there is a Kerdock code which is a 
non-linear OA(2 2 ™,2™,2,5). It has | as many points as the 
corresponding [2 m ,2m+ 1,6*] 2 BCH code flUtll- For any 
even m > 6, there is a Delsarte-Goethals code which is a non- 
linear OA(2 3m - 2 ,2"\2,7) H. It has \ as many points as the 
corresponding [2 m , 3m + 1 , 8*] 2 BCH code. 

The following result comes from thinning some cubature 
formulas of Stroud, some of whose points have a product 
structure. 

Theorem 3.1. Let n>6 and let 

{4m 2 2m - l <n<2 2m 
4m + 2 2 2m <n<2 2m + 2 m 
4m + 3 2 2 "'+2 m <«<2 2m+1 

Then the sphere S 1 " , R n with Gaussian measure, and the ball 
B n admit QLPI 5-cubature formulas with 2 k + 2n points. 

2 

Proof. The formulas S n :5-3, U„:5-2 and E' n :5-3 listed in 
Stroud 1 20, pp. 270,294,317] have 2" + 2n points with 2" of 
them lying on the vertices of a cube. These 2" points can be 
thinned to either the [2 2m+1 ,4m + 3,6*] 2 BCH code, or the 



Actually, Theorem O.ll is not quite optimal, because it uses 
a convenient set of good distance-6 linear codes and non- 
linear strength-5 orthogonal arrays rather than the best ones 
presently known. A complicated map of the best presently 
known linear codes over F 2 of length n < 256 is provided by 
the "best codes" functions in Magma |22]. Undoubtedly this 
map could be augmented by non-linear orthogonal arrays, but 
we know of no effort to do so. When n is a power of 2, Ker- 
dock and BCH codes are the best presently known choices. 

Victoir 12111 also established Theorem 13.11 (with BCH 
codes). If n = 2 m and Stroud's formulas for S' 1 ^ 1 is thinned 
using a BCH code, it then has equal weights and is therefore a 
5-design. Interestingly, in this case it has a transitive symme- 
try group and was previously found by Calderbank, Hardin, 
Rains, Shor, and Sloane |3[ . Similar constructions were found 
by Konig 1 12], by Sidelnikov fl9ll . and by Schechtman, inter- 
preting work of Hajela @]. 

We can obtain a good 3-cubature formula for the cube C„ by 
a straightforward application of Theorem ll.2l using the 2-point 
Gaussian quadrature formula for the interval [—1,1]. Thinning 
the product formula using a BCH code yields a 2 ;+1 -point 
formula when 2 ; '~ 1 < n < 2 J . When n — 2 J , or more generally 
whenever there is a Hadamard matrix of order n, the product 
formula can be thinned to the 2n vertices of a certain regular 
cross-polytope inside C„. A formula due to Stroud (C„:3-l 
12(1 p. 230]) also uses the vertices of a regular cross-polytope, 
but not the same one. 

We can obtain a 3-cubature formula with 0{n) points for 
A„_i with a similar construction. Using known Hadamard ma- 
trices, the formula has 3n + o{n) points; if the Hadamard con- 
jecture holds, it has between 3n — 1 and 3n + 5 points. First, 
the positive ray R>o with exponential weight has a equal- 
weight 2-quadrature formula with points at and 2. If we ap- 
ply Theorem l 1.1 I to this formula and a Hadamard matrix of or- 
der n, the result is a 2n-point formula F on K." which also has 
degree 2. However, if our interest is integration on A„_i, we 
need only consider homogeneous polynomials on M" . The 
formula F correctly integrates every degree 3 monomial other 
than x^. We can fix F for these monomials, without changing 
its sum for xfxj oxxiXjXk, by adding a point at (1,0,0,...) (and 
permutations) with weight 2. 

The projected formula on A„_i consists of these points and 
weights in barycentric coordinates: 



(«'«>•••' „>0i0>---)0)# 2 („+i)(„ +2 ) ■ 
(- - -) 4n 

Vn'n'"''n/ («+l)(n+2) 
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The subscript "S" denotes full symmetrization, as in Stroud's 
notation. The subscript "H" denotes symmetrization in the 
pattern of a Hadamard design. (See Section 0) This pro- 
duces a formula with 3« — 1 points provided that there exists 
a Hadamard matrix of order n. When there is none, we can 
use a Hadamard matrix of order £ > n. The formula on A^_i 
with 3£ — 1 points can be projected onto A„_ i , as in the proof 
of Theorem II. 31 We can take I = n + o(n) by letting £ — 1 
be the first prime after n which is 3 mod 4. If the Hadamard 
conjecture holds, we can take £ = 4\n/4\. 

Stroud asked for a practical, PI 5-cubature formula for Cioo- 
Following Theorem 11.31 we can find one by thinning the 
product formula coming from the 4-point Chebyshev quadra- 
ture on [—1,1]. This product formula is the convolution of 
200 pairs of points, so we can thin it using the Kerdock 
OA(2 16 ,2 8 ,2,5), projected to 200 dimensions. The cubature 
formula therefore has 2 16 = 65536 points, which would have 
been fairly practical even in 1971 when Stroud asked the ques- 
tion. (The Kerdock code used here was discovered shortly af- 
terward Ql but the BCH codes was known in 1959 lH.IToll.') 

Victoir 12111 found another thinning of the same Chebyshev 
product formula with 4 12 = 16777216 points, which the au- 
thor tied in the first version of this paper. 

Note that the Chebyshev-Kerdock 5-cubature formula for 
Cioo is overdetermined. The threshold of exact determina- 
tion for centrally symmetric 5-cubature formulas on Cioo is 
87651 points. Meanwhile the centrally symmetric Tchakaloff 
bound is 8852652 points, while the Stroud lower bound is 
5050 points. 

Finally Schiirer fl7ll compared the numerical accuracy of 
various cubature and quasi-Monte-Carlo methods for the inte- 
gration of various test functions defined on C„ with 2 < n < 
100. He assumed a more modern limit of 2 25 evaluations of 
the integrand. For much of this test regime we can suggest the 
following cubature formulas: Start with the power of the con- 
volutional 7-quadrature formula fl3ll for [—1,1], whose points 
are approximately at 

±.500128 ± .243941 ± .153942. 

Then thin the n-fold product power of this formula using a 
Delsarte-Goethals code. The result is an EI 7-cubature for- 
mula with at most 2 23 points up to dimension [256/3J = 85. 

4. OTHER COMMENTS 

Victoir l2lll proposes thinning symmetric cubature formu- 
las rather than product or convolution formulas. The enabling 
result of symmetric cubature formulas is Sobolev's theorem: 



If a linear action of a finite group G preserves /J., then a cu- 
bature formula consisting of orbits of G need only be checked 
for G-invariant polynomials. Victoir extends Sobolev's the- 
orem with a G-invariant generalization of Tchakaloff 's theo- 
rem: A PI cubature formula only needs as many orbits as the 
dimension of R[x]< ( , the space of G-invariant polynomials of 
degree at most t. One important special case is when G is 
the 2-element central symmetry group. If ju is a measure on 
W with central symmetry and t is odd, the bound from this 
version of Tchakaloff 's theorem is 0(n'~ v ) points. 

Even if a cubature formula F uses very few orbits of G, 
some of the orbits might be very large. Victoir proposes thin- 
ning each large orbit separately. He notes that this can be 
done using linear programming, among other methods; linear 
programming on a set of G-orbits should be much easier than 
general numerical methods to find positive cubature formulas 
for /i. If G = (Z/2)" is the group of independent sign changes 
of all n coordinates, then an orbit of G is a Cartesian power 
can be identified with F| for some k <n. In this case Victoir 
found the constructions of Theorem l 1.11 and Theorem ll.2l (In 
the case of 5-cubature on C„, he found a special construction 
with 0(n 3 ) points with elements of both Theorem l 1.1 1 and the 
n-cube case of Theorem ll.3P 

If G is the group of coordinate permutations, then an orbit 
whose points have two distinct coordinates can be identified 
with the set of A:-subsets of an n-set. A geometric f-design T 
within this orbit is also a traditional combinatorial f-design, 
or an (n,k,t) — X design. Namely, T is a collection of blocks 
of size A: in a set of n such that each f-subset is contained 
in exactly X blocks. In particular, an (n, |,3) — | design is 
called a Hadamard design, because it comes from the rows of 
a Hadamard matrix. 

These constructions motivate the notion of a weighted or- 
thogonal array. We define it as a finite set A and a measure 
jj. on A" that projects to uniform measure on each A 1 with 
|/| < t . More generally, jj, might project to a" for some ref- 
erence measure a on S. Such arrays could improve of The- 
orem the factor formulas would not need to have equal 
weights. 

Finally, cubature formulas coming from Theorem l 1 . 1 I could 
be viewed as quasi-Monte-Carlo methods. They are simi- 
lar to some constructions of (f ,m,i)-nets, which are quasi- 
Monte-Carlo methods first defined and largely developed by 
Niederreiter 11511 . Nonetheless PI cubature formulas and 
discrepancy-based quasi-Monte-Carlo methods are thought to 
have complementary advantages fl7ll . We believe that the im- 
proved asymptotics presented here could change the standing 
of cubature among numerical methods for integration. 
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